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Abstract 


This report Is concerned with the derivation of the equations of 
motion for the Spacecraft Control Laboratory Experiment (SCOLE). ."or 
future reference, the equations of motion of a similar structure 
orbiting the earth are also derived. The structure is assumed to 
undergo large rigid-body maneuvers and small elastic deformations. A 
perturbation approach is proposed whereby the quantities defining the 
rigid-body maneuver are assumed to be relatively large, with the elastic 
deformations and deviations from the rigid-body maneuver being 
relatively small. The perturbation equations have the form of linear 
equations with time-dependent coefficients. An active control technique 
can then be formulated to permit maneuvering of the spacecraft and 
simultaneously suppressing the elastic vibration. 
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1, Introduction 


We shall first derive Lagrange's equations of motion for the 
spacecraft of Fig. 1 regarding the structure as orbiting about the earth 
and then modify these equations so as to describe the laboratory 
experiment. In the derivation, the shuttle is treated as a rigid body 
and the beam and antenna as flexible, distributed parameter systems. 

The equations of motion can be further modified for the case of a rigid 
antenna. 

The equations describing a maneuver of a rigid space structure 
consist of nonlinear ordinary differential equations. On the other 
hand, the equations describing the small elastic displacements of a 
flexible structure, relative to the rigid-body maneuver are linear 
partial differential equations. Hence, the complete equations of motion 
describing a flexible bo u y during a maneuver represent a set of 
nonlinear hybrid differential equations. 

Hybrid systems possess an i nf ^ L number of degrees of freedom. 

In practice, however, it is necessary to reduce the number of degrees of 
freedom to a finite one, which implies spatial discretization and 
truncation. Substructure synthesis often proves useful as a method of 
discretization and truncation, particularly in the case of dL— ibuted 
substructures. Even in the case of discrete substructures, et of 

linearly independent vectors can be used as admissible vectors to reduce 
the number of equations of motion. 

In this report, we propose a perturbation technique whereby the 
flexible spacecraft maneuver is assumed to consist of a combination of a 
rigid-body maneuver and small motions including rigid-body deviations 
from the rigid-body maneuver and elastic vibrations. Regarding the 



rigid-body maneuver as known, the perturbation equations for the 
vibration control reduce to a set of linear ordinary differential 
equations with known time-varying coefficients. 

2. Equations of Motion of the Spacecraft 
It is convenient to refer the motion of the spacecraft to a given 
reference frame XQy 0 z 0 , where the frame can be regarded as being 
embedded In the rigid shuttle. The reference frame has six degrees of 
freedom, three rigid-body translations and three rigid-body rotations. 

We propose to derive the equations of motion by means of the 
Lagrangian approach. To this end, we must first obtain expressions for 
the kinetic energy, the potential energy and the virtual work. 
Considering Fig. 1 and denoting the position of the origin 0 of the 
frame Xgy 0 z 0 b Y the vector ? and the position of a point S in the 
shuttle relative to 0 by r, the position of S relative to the inertial 
frame XYZ is R s = R + r. Moreover, denoting by a the vector from 0 to a 
nominal point A on the appendage and by u the elastic displacement 

vector of the point, the position of A In the displaced configuration 

is R. = R + a + u. It must be noted that the vectors r, a and u are 

" n ****** r» W 

likely to be measured relative to axes XqYqZq. In view of the above, 
the velocity of a point S in the shuttle is 
* * 

Rj = R + id x r (1) 

where R is the translational velocity and w is the angular velocity of 
the frame XoygZg with respect to the inertial frame. Similarly, the 
velocity of a point A in the appendage is 

R. = R + ui x (a + u) + u (2) 
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where u Is the elastic velocity of the point relative to the x 0 yQZg 
frame. Hence, the kinetic energy of the spacecraft Is 


T *|J I8$IS + R I 5a I Z dra A 

m s » A 

■ I J I? + “ « cl Zdn, s + 1 -f IB + s * (J + u) + Bl 2dm « 
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= I m|R| 2 + j w T I 0 u + R*(w x S Q ) + i J | u I' 
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+ R-[J u dm^ + u x J u dm^] + / u*(w x a)dm^ 
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+ J (w x a) • (w x u)dm^ "2 J | u> x u | dm. 
m, m* 
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+ J u * (u> x u) dm. 

m fC 


where 


( 3 ) 


s 0 - ; r dm s + J a dm. (4) 

"S ”A 

and m^ and m are the masses of the shuttle, appendage and entire 
spacecraft, respectively. Also, Iq is the total mass moment of Inertia 

n 

matrix of the undeformed structure about point 0. Note that |x| 
denotes the inner product x*x. 

The potential energy is due to the combined effects of gravity and 
strain energy. Assuming that the origin of the inertial coordinate 
system coincides with the center of the gravitational field, the 
gravitational potential can be expressed as 



cr* 

where m e Is the mass of the earth and G Is the gravitational constant. 

The strain energy can be expressed as an energy Inner product 
symbolized by [ , j (Ref. 1). The total potential energy then becomes 

' V = ^[u,u] 4- V g (6) 

The virtual work Is due to external forces, including control 
forces. Denoting by f s the force vector per unit volume of the shuttle 
and by the force vector per unit volume of the appendage, we can 
write the virtual work as 

*" * L fs' s Bs d0 s + T Ea-Va ‘ 7 > 

U S U A 

where D<* and are the domains of the shuttle and appendage, 
respectively. 

Before deriving the equations of motion, we consider certain 
simplifying assumptions. To this end, we estimate the maximum possible 
angular velocities by Ignoring the effects of the appendage and 
examining the prescribed maneuvers of the shuttle alone. For the 20° 
maneuver about the Xq axis, applying maximum torque, the maximum angular 
velocity is approximately .06 rad/s. For the 90° maneuver about the z 0 
axis the maximum angular velocity is .047 rad/s. If the elastic 
displacements are small, then the last two terms of Eq. (3) are of 
higher order and can be neglected. The third to last term of Eq. (3) 
will be retained, despite leading to nonlinear terms, because the 
magnitude of the factor multiplying the independent variables u and 
to tends to offset the smallness of the independent variables. Next, we 
express the elastic displacements in the form of linear combinations of 
admissible functions, or 
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u - *g (8) 

where t Is a matrix of space-dependent admissible functions and g is a 
vector of time-dependent generalized coordinates. Introducing Eq. (8) 
Into Eq. (3) and neglecting the last two terms in Eq. (3), the kinetic 
energy takes the matrix form 

T 5 i n,R T R ♦ i S T I 0 « + R T C T S oS * \ g\g 

+ R^C^$g + R^c\J$g + gWu + a^w^o dm.g (9) 

~ \ 

where 

, T T 

* = J « dm., $ = J* ra dm A {10a, b) 

m A m A 

and 

M a = J o T « dm A (10c) 

m ft 

is the mass matrix of the appendage. The symbol C represents a rotation 
matrix from the Inertial frame to the x 0 ygz 0 frame and Its elements are 
nonlinear functions of a set of Euler angles a. The tilde over a given 
vector such as v denotes a skew symmetric matrix of the form 



Recognizing that the magnitude of R is large and u is small in 
comparison with the other vectors in Eq. (5) and ignoring terms of order 
higher than three, a binomial expansion permits us to write 

\ s - GHeMSl ' 1 - 5-(? 0 + I u dm.)|R|- 3 l (12) 

"A 
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Introducing Eq. (8) Into Eq. (6) and considering Eq. (12), the potential 
energy can be written in the matrix form 

, T Gm_m Gm T t _ 

v ■ ? 9 k a9 " JET + - C (? o + *9) < 13) 

where 


K a = l®, ®] (14) 

Is the stiffness matrix of the appendage. The virtual work can be shown 
to have the expression 

6W = F T CsR + M^Sa + Q^sq (15) 

where 


E - J„ £s dD s + Ja dD A 

U S U A 

5 = L dD S + L *\ dD A + T. f A 4 dD A9 ( 16 > 

U S U A U A 

9 ■ t *\ dD A 

U A 

are generalized force vectors In terms of components about xq,y 0 and z Q . 

Without loss of generality, we let point 0 correspond to the center 
of mass of the spacecraft ‘in its undeformed state, so that the vector 
Sq Is zero. Then, Lagrange's equations of motion can be written in the 
symbolic form 


a! (4) ♦ S! - c t f 


'aft 


(17a) 


A (ill 

dt 'v J 
3a 




(17b) 


_d f aT\ IT 
dt aq " 8 9 



(17c) 
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t 


so that, considering Eqs. (9), (13) and (15), the equations of motion 


for the spacecraft in orbit are 


S „T— •• „ , ”5 cT »3 R(3R T C T fg) T 
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t Or _ 

IqU + w Iqw + ( CR ) + w 5g + J(u>)g + J(w)g 

x Gro _ 

+ 5 J(w)g + ~~3 [CRjtg 3 M 

M Grci 

H ft g + <t»"^CR + 5^w + « A g + J c^w^a dro^w + — ^ $^CR = Q 

m„ R I 


(18a) 


(18b) 

(18c) 


where 


J(w) a J (aS + [au])« dm A (18d) 

m A 

Higher-order terms have been neglected In Eq. (18b), consistent with the 
preceding discussion concerning the magnitude of the maneuver angular 
velocities. However, as before, nonlinear terms with large 
coefficients, such as Iq and am^, have been retained. The position 
vector R, Its time derivatives, and the Euler angles vector a have been 

f* ~ 

considered to be of arbitrary magnitude, with the result that many other 
nonlinear terms have been retained in Eqs. (18). 

3. Equations of Motion for the Laboratory Experiment 
In the laboratory experiment, the spacecraft Is not actually free 
in space, but suspended from the celling by means of a cable or a 
beam. The following analysis applies to either case. The support Is 
likely to affect the dynamic characteristics of the system. Hence, In 
the sequel, the support is added to the free model in the form of an 
elastic member. 
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Considering Fig, 2, the position vector for an arbitrary point C on 
the cable Is R r = c + w, where c is a position vector and w Is the 
elastic displacement of the cable, both of which are measured with 
respect to the inertial frame. The position vector for the point 0 Is 

R = c B + w Q + e (l f j) 

where the subscript B denotes evaluation at the point B and e Is the 
vector from point B (ball joint) to the point 0 fixed on the "shuttle”, 
measured with respect to the XqYqZq frame. The velocity vector of an 
arbitrary point C on the cable Is then 
« 

Be ■ S (20) 

and the velocity of point 0 Is 

t 

R = W B + u X e (21) 


The kinematics for the shuttle body and appendage remain the same as for 
the unrestrained spacecraft In space. Hence, the kinetic energy for the 
entire structure is 


T = ij |R c | 2 dm c + |j |R s | 2 dm s +ij |R A | 2 dm A 
m c m $ m A 

■if | w 1 2 dm c + A J |w B + oj x (e + r) | 2 dm- 

m C ~ m S ~ 

+ 7T / |w B + tox(e + a + u)+ u| dm A 
m A 

= \ l |w| 2 dm c +|m|w B | 2 

m c 

+ W B -(w x S B ) + (w X e)* (at X Sg) + ~ f |u| 2 dm A 

m A 

+ w B - J u dm A + (u x e) -J u dm A 
ffl a m a 
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( 22 ) 


<* 


+ j‘ u* (u « a) dm, + / (aj X a)'(w X u) dm A 
\ ~ " m A 

where 

r B ■ J 0 + m5Tg * 5 b “ §0 + < 23 > 

The elastic displacement vectors u and w have been assumed to be 
small. Also, the angular velocity was assumed to be, at most, of the 
same order of magnitude as In the case of unrestrained spacecraft, so 
that terms of higher order can once again be neglected. 

The expression for the virtual work Is given by Eq. (7), but the 
potential energy must be modified. The acceleration of gravity will be 
assumed to be constant and can be expressed as g = - gu^ so that the 
gravitational potential Is 



where R Is defined by Eq. (19) and u. Is a unit vector In the Z 

N ** Im 

direction. The elastic potential energy for the system can be expressed 
as 

V, = ~ lu,u] + A [w,w] (25) 

L C ~ ~ C ~ ~ 

where [w,w] is the energy inner product for the cable, which includes a 
stiffening effect due to the weight of the "spacecraft". 8ecause a 
cable has little inherent bending stiffness, this effect can be 
significant. 

As with the appendage In the preceding section, the elastic 
displacement of the cable can be approximated by a linear combination of 
admissible functions, or 
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w ■ i»n 

M * W 


( 26 ) 


1 . 


e 


where t|» Is a matrix of space-dependent admissible functions and n is a 
vector of time dependent generalized coordinates.' Introducing Eqs. (§) 
and (26) Into Eq. (22), the kinetic energy takes the matrix form 



where 

M c = / dm c + nukgijig (27a) 

m c 

is the mass matrix of the cable, including the mass of the structure 
lumped at the end of the cable. ^ denotes the matrix ip evaluated at 
b. Introducing Eqs. (8) and (26) Into Eqs. (2d) and (25) the potential 
energy can be written In the matrix form 

V = + sjcg + mn T ^g + g Vc 3 + 1 gT^ + 1 ^ (ZB) 

where 


? = J ^ dm c (29a) 

m C 
and 

K c = (29b) 

Is the stiffness matrix of the support. Considering Eqs. (19) and (25), 
the virtual work can be expressed as 

sW = F^CiMn + M*^6a + Q^dq (30) 

M Q fV ~ Mf «V 
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I 

where all the terms have been defined previously. The effect of the 
friction of the ball joint can be assumed In the form of an external 
torque. Hence, we let M = M_ + M* In Eq. (30), where Is a vector of 
control moments and M f Is a vector of frictional moments caused by the 
ball joint. 

Without loss of generality, we assume that point 0 Is the center of 
mass of the spacecraft, so that Sg Is zero. For the laboratory 
experiment to be successful the vector e must be very small, so that the 
center of mass 0 tends to coincide with point B. This Is difficult to 
accomplish in practice due to the friction of the ball joint and the 
difficulties of dynamically balancing a large flexible structure. 
Realistically, we can expect e to be small but nonzero. Hence, for 
simulation purposes, we choose e to be small, so that the kinetic energy 
of Eq. (27) can be approximated as follows: 


T 5 i 2 Tm c3 + \ + m D T 1 J B cT ^ + \ 9\9 


T T T T T T T 

r n ijjgC ‘sg + cjOu + to J* a w $ dm^ g 

m A 


(31) 


Lagrange's equations remain in the symbolic form of Eqs. (17), with 
the exception of Eq. (17a) which must be replaced by 


d ( 3 l\ , aV . . T r Tr 


(32) 


This results in the following equations of motion for the laboratory 
experiment: 


JnTrT, 


JrTa- , Jr Tt: 


J . Js 


M c n + mipgC'iu'ew + mipgC'eu + nigC'^g + (mtjjg + «j»')g + K c n = 4 »gC T F 

(33a) 


Iq(d + + 5>g + u 1 5g + J(w)q + J(u))g + uj*J(uj)q 


~T~ 


~T • 
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(33b) 


fro 

•p H 

+ e Cipgnm + {Cg ) (em + ®g) * H 

Mg + 5 u + © C^gn + J © w a dm. w + © Cg + Kg - Q (33c) 

; m A ~ 

Higher-order terms have been neglected and some nonlinear terms retained 

in Eq. (33b), consistent with the case of the spacecraft in space. In 

this case, the Euler angles vector a is assumed to be of arbitrary 

magnitude, which Is responsible for many nonlinear terms in Eqs. (33). 

4. Simulation and Contr ol 

The nonlinear equations of motion take the same basic form for the 
orbiting spacecraft as for the laboratory experiment, as can be seen 
above. Hence, the solution techniques suggested here apply for both 
situations. 

Consider a first-order perturbation in the vectors R and u. 

5 “ 5o + h* ~ = So + ~i ( 3 4) 

where the first-order terms and are small compared to the zero- 
order terms R Q , to Q . Introducting Eqs. (34) Into the nonlinear equations 
of motion, Eqs. (18) or (33), we obtain zero-order and first-order 
equations. The zero-order solutions R 0 and oj q can be obtained from an 
open-loop, rigid body maneuver strategy. These solutions can then be 
inserted into the first-order equations yielding linear equations with 
known time-varying coefficients. Simulation requires numerical 
Integration. The equations can be put into their most compact form by 
approximating the first-order motion by means of a linear combination of 
eigenvectors corresponding to the lower frequencies. A control 
technique suppressing elastic vibration as well as rigid-body motion 
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that deviates from the desired maneuver, can also be formulated 
considering the equations In this compact form. 

5. Conclusions 

i 

The equations of motion for the structure both In orbit about the 
earth and In the laboratory are nonlinear, even when the elastic 
deformations are small. The nonlinear terms result from the large 
rigid-body maneuver. For instance, centrifugal and Coriolis terms 
affecting the elastic displacements can be significant in a minimum-time 
maneuver, so that they must be retained In the equations. Through a 
perturbation approach, the nonlinear equations of motion can be 
transformed Into a set of equations governing the rigid-body motions and 
a set of time-varying, linear equations governing small deviations from 
the prescribed rigid-body maneuver and elastic motion. Future work will 
include applying this perturbation approach to the above equations of 
motion for use in both simulation and control of the structure. 
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